Electron-phonon interaction in the t-J model 
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We derive a, t-J model with electron-phonon coupling from the three-band model, considering 
modulation of both hopping and Coulomb integrals by phonons. While the modulation of the 
hopping integrals dominates, the modulation of the Coulomb integrals cannot be neglected. The 
model explains the experimentally observed anomalous softening of the half-breathing mode upon 
doping and a weaker softening of the breathing mode. It is shown that other phonons are not 
strongly influenced, and, in particular, the coupling to a buckling mode is not strong in this model. 



There has been a strong interest in electron-phonon 
coupling for high- Tc cuprates after the discovery of strong 
coupling to a mode at 70 meV in many cuprates |]| . The 
coupling was ascribed to a half-breathing phonon along 
the (1,0,0) direction, i.e., an in-plane bond-stretching 
mode, with an electron-phonon coupling A ^ 1. This 
phonon shows an anomalous softening when the cuprates 
are doped, in particular, towards the zone boundary 
I2, y, y, I3 ■ While the softening of other phonons upon 
doping can be understood as a screening of the ions, the 
softening of the half-breathing mode cannot be described 
in a shell model with conventional parameters This 
suggests a substantial electron-phonon coupling. The 
phonon has an appreciable broadening 0, also support- 
ing a substantial coupling. Local-density approximation 
(LDA) calculations, however, show a weak coupling to 
this phonon, with A at the zone boundary being ~ 0.01 
Anomalous behavior of bond-stretching modes has 
also been observed in other compounds 0, 0, E| ■ It is 
then interesting to study the electron-phonon interaction, 
taking into account the strong electron-electron interac- 
tion 12 0. This has been done in the t-J model 
by von Szczepanski and Becker by KhaliuUin and 
Horsch and by Ishihara, Nagaosa and coworkers . 

The t-J model is derived from a three-band model of 
a Cu02 layer including a Cu dxi^yi and two O p 
orbitals. To obtain the electron-phonon coupling, the 
atoms are displaced. This leads to a change of the hop- 
ping integral tpd and the charge transfer energy Sp be- 
tween Cu and O atoms, where the change of £p is assumed 
to be due to a change of the Coulomb integral Upd be- 
tween nearest neighbor Cu and O atoms. Transforming 
to a, t-J model, this leads to changes of both the on-site 
energy Eg of the Zhang-Rice singlet and the hopping be- 
tween the sites (off-site coupling) . Two groups I4, il5(| 
assumed the variation of tpd to dominate, while a third 
group 16] focused on Sp. The first two groups assumed 
the off-site coupling in the t-J model to be negligible, 
while the third group emphasized this coupling. 

The singlet energy in the t-J model is large, \Es\ ^ 5 
eV. For a rigid lattice and fixed doping, this energy enters 
only as an uninteresting constant. Because of the strong 
distance dependence of tpd^ however, we may expect a 
strong electron-phonon coupling from this term. 



We study modulations of tpd and ep by phonons. 
The tpd modulation dominates, but destructive inter- 
ference between the two effects is not negligible. The 
phonons couple mainly to on-site terms. Phonon spec- 
tral functions are obtained from exact diagonalization. 
The model explains the anomalous softening of the half- 
breathing mode and correctly gives a weaker softening of 
the breathing mode. Other phonons are only softened 
weakly. Experiments suggest a strong coupling of prefer- 
entially the anti-nodal electronic state to a phonon at 40 
meV, perhaps a Big buckling mode 0| . We find a weak 
coupling to the Big mode in the model studied here. 

The derivation of a model with phonons proceeds as 
for the normal t-J model but with displaced atoms. 
With the atoms in their equilibrium positions, a Cu atom 
couples only to a given linear combination of p-orbitals on 
the four O neighbors 0|. With displaced atoms there is 
also coupling to other combinations. This coupling, how- 
ever, enters to second order in the displacement and is 
neglected here. Following Zhang and Rice J^j, we work 
to lowest order in tpd, although this is not very accurate 
for realistic parameters. We neglect certain terms where 
a hole hops i k j {k ^ i, k ^ j). We also neglect 
the electron-phonon interaction via superexchange inter- 
action, since this coupling occurs only via motion of Cu 
atoms and with a small prefactor ■ We then obtain a 
t-J model with the linear electron-phonon coupling 
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Here cj^ creates a d-hole on site i if this site previously 
had no hole, and 5jj ^ creates a phonon with wave vector 
q, index v and frequency uji, (q) . The system has N sites 
with coordinates R^. (5i.j-x = 1 if site j is one site to the 
right of i in the x-direction. The modulation of tpd gives 
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with an analogous form for . Here s; = 
ci = cos{qia/2), l^x,y\X = X^/^k = 7o-7i = 0.96, 
ri = 4(72 - 70), T2 = 871 - 70 - 572 - 2711 and 
rs = 871 - 70 - 72 - 6711, with 70 = 7(0) = 1.285, 

71 = 7(ax) = 0.327, 711 = 7(ax + ay) = 0.209 and 

72 = 7(2ax) = 0.164, where 7(R) = ^kc''"-^/^ with 

/3k = l/y's^ + and a the lattice parameter. Mcu and 

Mo are the masses of the Cu and O atoms, respectively, 
Eq^ is the x-component of the polarization vector for 
one of the O atoms (O^;) with an analogous meaning of 
Cq^. U is the hole-hole repulsion on Cu. von Szczepanski 
and Becker calculated gijiq, v) for the q — (1, l)7r/a 
breathing mode and found (q, v) = Q (since ef,^ = 
for this value of q), as above, and a slightly different re- 
sult for A(q, v) due to slightly different approximations. 
Considering a variation of the charge transfer energy, 
assuming that it is driven by a nearest neighbor Cu-0 
Coulomb interaction Upd, we find 
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Here a = J2 [3^^^ e'""^-^ / N = -0.14. 

To estimate the relative magnitude of these terms, we 
use [12, 20] tpd = 1.2 eV, £p = 3 eV, f7 = 10 eV, Upd = 1 
eV, and a — 3.8 A. We assume the distance dependence 
tpd ~ with n = 3.5, based on LDA calculations [2lj |. 
and Upd ^ r^^. This gives the hopping integral -0.47 
eV in the t-J model, close to values commonly used [2^ . 
For q = (g,,0), a; = 0.07 eV, eg, 
polarization vectors zero, we obtain (in eV) 



1 and all other 
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FIG. 1: Phonon dispersion in the (1,0) and (1,1) directions. 
Experimental results (dotted line) for a; = and x — 0.15 are 
shown. Theoretical results (full curve) for x — 0.125 show 
the calculated softening from the experimental a; = results. 
The average over boundary conditions is shown and the bars 
show the spread due to different boundary conditions. There 
is a strong softening in the (1,0) direction, while the softening 
in the (1,1) direction is weaker. 



due to modulation of Upd- This shows that the modula- 
tion of tpd dominates over that of Upd- This is mainly due 
to a stronger power dependence (n = 3.5) for tpd than 
Upd but also due to a partial cancellation between the 
terms proportional to l/e^ and 1/{U — e^)^ in the con- 
tribution from dUpd/dr. The contribution from dUpd/dr 
alone is very small. There is a destructive interference 
between contributions from dtpd/dr and dUpd/dr, how- 
ever, which reduces the phonon softening by about 30%. 
The prefactors of the off-diagonal terms are small [Eqs. 
Q and ((SJ], and the diagonal term dominates. We have 
neglected quadratic terms in the phonon displacement in 
ljT]l , although some terms may give non- negligible contri- 
butions to the doping dependence of phonon energies. 

The model describes the softening of phonons due 
to holes in the doped system, but it does not include 
other interactions present in both the doped and undoped 
systems. These interactions are described by a two-spring 
model, fitted to the phonon frequencies in the (1,0) and 
(1,1) directions of the Mnrfopec? system. This spring model 
provides the eigenvectors e in Eqs. (PJ-®. 

To study the t-J model with phonons, we use exact di- 
agonalization |2^ , including all possible electronic states 
for a finite cluster of size M x N. To obtain a finite 
Hilbert space, we allow states containing a maximum of 
K{= 5) phonons '2^1 , which is sufhcient for convergence. 
We calculate the phonon spectral function {u; > 0) 

i?..(c.) = -ilm(0|V^;;-^^^^^^<jO), (9) 



where |0) is the ground-state with the energy Eq, i/iq^ — 
bqu + b^-qy and 77 is infinitesimal. Since the clusters that 
can be treated are too small to give a quasicontinuous 
spectrum, we use the center of gravity of the phonon 
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FIG. 2: Phonon dispersion in the (1,0) direction for clusters 
of different sizes. The corresponding doping is indicated. The 
figure shows that the softening increases with doping. 
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spectrum to define the renormalized phonon frequency. 
This provides the phonon softening in relation to the ex- 
perimental phonon frequencies of the undoped system. 
We have used the parameters above and J/t = 0.3. 

Figure ^ shows results in a 4 x 4 t-J model for the 
(1,0) and (1,1) directions. The dotted curves show exper- 
imental results 0, and the full curve shows the soften- 
ing due to the electron-phonon interaction in the doped 
system. Because of the small cluster sizes, the results 
depend on the boundary conditions. We have used pe- 
riodic, antiperiodic, and mixed boundary conditions, ap- 
plying periodic boundary conditions in one direction and 
antiperiodic in the other. Figure ^ shows the average 
and the bars the spread of the results. Although the re- 
sults are fairly sensitive to the boundary conditions, the 
trends are clear. Doping leads to a pronounced softening 
of the half-breathing mode along the (1,0) direction for 
|q| ""/(So), while the softening is weaker in the (1,1) 
direction, in agreement with experiment . Similar con- 
clusions have been obtained using analytical treatments 

The prefactor of the dominating diagonal term [Eqs. 
(O and lOl] in the coupling for the half-breathing mode is 

sin((72;a/2), suggesting a softening siv?{qxa/2). This 
behavior is essentially found in the calculations, although 
the softening is stronger for qx — 7r/(2a) than would be 
expected from this argument. The prefactor, however, is 
a factor \/2 larger for the breathing mode [q — 7r/a(l, 1)] 
than for the half-breathing mode [q = 7r/a(l, 0)], suggest- 
ing twice as large a softening for the breathing mode. Ac- 
tually, the softening is larger for the half-breathing mode, 
since it couples to excitations at lower energies. 

Figure 13 shows results in the (1,0) direction for differ- 
ent cluster sizes. Since all clusters have two holes, the 
change of cluster size changes the doping. The softening 
increases with doping x, e.g., at the zone boundary, the 
softening is 10% for x = 0.1 and 19% for x = 0.167. This 
is also observed experimentally, e.g, for La2-2;Sr2,Cu04 
the softening is about 10% for a; = 0.10 Q and about 14% 
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FIG. 3: Phonon dispersion in the (1,0,0) direction for a shell 
model. The dashed curves show results for the undoped sys- 
tem and the full curves results with the extra 0-0 spring 
describing the O atom coupling to the Zhang-Rice singlets. 
The arrows indicate the strong softening of the half-breathing 
mode, while other modes are not changed very much. 



for X — 0.15 3J. For the breathing mode [q — 7r/a(l, 1)], 
we obtain the softening 8% [x = 0.125), compared with 
the experimental results 3% {x = 0.1) and 8% (x — 0.15) 
i. 

The calculations show that there is a strong coupling 
to the half-breathing mode, where two O atoms move to- 
wards the Cu atom in between. At the zone boundary, 
the Cu atoms do not move. Towards the zone center, 
however, there is a substantial movement of Cu. Normal- 
ization of the polarization vectors then leads to a reduc- 
tion of the O polarization vectors. Completeness requires 
that the missing O weight is transferred to other modes. 
In a two-spring model the weight goes to an acoustic 
mode. Because of its small frequency, this mode is then 
softened far more [almost 50% for q = tt/ (2a)(l, 0)] than 
observed experimentally. 

To address this, we have used a more realistic shell 
model 2] for obtaining eigenvectors. This model gives al- 
most exactly the same eigenvectors for the half-breathing 
mode as our two-spring model. The "missing" O weight 
towards the zone center, however, is distributed over sev- 
eral modes, and the softening of a given mode is weaker. 
For instance, the longitudinal acoustic q = 7r/(2a)(l,0) 
phonon is softened by about 25%. Although smaller than 
in the two-spring model, the softening is still too large. It 
is, however, further reduced by the repulsion from lower- 
lying modes of the same symmetry. 

To study this, we have modified the shell model 2] to 
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take the electron-phonon interaction into account. The 
movement of two O atoms towards a Cu atom leads to 
a lowering of the Zhang- Rice singlet energy. The system 
can take advantage of this by transferring a singlet to 
such a site. This is approximately described by intro- 
ducing a spring with a negative spring constant, k = — 3 
eV/A^, between two O atoms on opposite sides of a Cu 
atom. A similar term was used to describe La2-£!;Sra;Ni04 
0, . The present work gives justification for such a 
spring. Figure 13 compares results of the shell model with 
and without the additional spring. Apart from the half- 
breathing mode, no mode is strongly softened by the new 
spring 26|. The t- J model thus correctly softens the half- 
breathing mode, without introducing unphysical soften- 
ing of other modes. 

An LDA calculation for the frequency of the half- 
breathing mode in (doped) YBa2Cu307 found good 
agreement with experiment although the very small 
calculated electron-phonon coupling would suggest a 
weak doping-dependence. Since LDA cannot describe the 
insulating undoped system, the implications for the dop- 
ing induced softening are unclear. 

The half-breathing and (in particular) breathing 
modes have unfavorable q dependences for d-wave pair- 
ing. The coupling to these modes is not expected to 
explain superconductivity in Eliashberg-like theories. 

Recent photoemission experiments suggest strong cou- 
pling to a phonon at 40 meV, perhaps a Big buckling op- 
tical phonon, i.e., a phonon where the in-plane O atoms 
have c-axis out-of-phase displacements For symme- 
try reasons, this mode couples only to second order for 
a single flat Cu02 plane, due to the hopping between 
the Cu d^2_y2 and O Pz orbitals. Such terms, however, 
have a small prefactor, and the coupling constant for the 
on-site term is more than an order of magnitude smaller 
than for the half-breathing mode [l^. The planes may 
have local static bucklings, allowing coupling to linear 
order. Assuming a 0.2 A (6°) buckling, we find that the 
coupling constant is still an order of magnitude smaller 
for the buckling mode. It seems that one would have to 
go beyond a single layer t-J model to obtain a strong 
coupling to the buckling mode. 

We have found that the distance dependence of tpd is 
substantially more important than that of Upd for the 
electron-phonon interaction. Nevertheless, the interfer- 
ence effects cannot be neglected. The t-J model with 
phonons describes the strong renormalization of the half- 
breathing mode, a weaker renormalization of the breath- 
ing mode and no anomalies in other modes. 
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